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Abstract. If dark matter is unstable and the mass is within GeV-TeV regime, its decays pro¬ 
duce high-energy photons that give contribution to the extragalactic gamma-ray background 
(EGRB). We constrain dark matter decay by analyzing the 50-month EGRB data measured 
with Fermi satellite, for different decay channels motivated with several supersymmetric sce¬ 
narios featuring R-parity violation. We adopt the latest astrophysical models for various 
source classes such as active galactic nuclei and star-forming galaxies, and take associated 
uncertainties properly into account. The lower limits for the lifetime are very stringent for a 
wide range of dark matter mass, excluding the lifetime shorter than 10 28 s for mass between 
a few hundred GeV and ~1 TeV, e.g., for bb decay channel. Furthermore, most dark matter 
models that explain the anomalous positron excess are also excluded. These constraints are 
robust, being little dependent on astrophysical uncertainties, unlike other probes such as 
Galactic positrons or anti-protons. 
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1 Introduction 

Recent cosmological measurements of the cosmic microwave background (CMB) [1, 2] con¬ 
clusively show that about 84% of matter in the Universe is made of dark matter. Evidently 
these dark matter particles are electromagnetically neutral and also stable compared to cos¬ 
mological time scale. However, it does not necessarily mean that they have infinite lifetime. 
In most cases, dark matter is simply assumed to be stable, e.g., by introducing Z 2 symmetry 
in the Lagrangian. One of the most popular examples of such symmetries is R-parity in 
supersymmetric models, which is often introduced to make protons stable. It also stabilizes 
the lightest superparticle (LSP) to be a candidate for dark matter. However, whether dark 
matter is absolutely stable or not, and hence whether such symmetry as R-parity exists or 
is broken, should be tested in light of modem data. 

Dark matter decays produce cosmic ray particles, and details of the final state de¬ 
pends on its interaction with standard-model particles. In supersymmetry, R-parity viola¬ 
tion (RPV) allows the LSP to decay to standard-model particles, and this is indeed favorable 
for some supersymmetric scenarios from cosmology arguments. For instance, if gravitino, 
a superpartner of the graviton, is the LSP, the next-LSP may decay during big-bang nu¬ 
cleosynthesis, which destroys light elements. In RPV scenario, however, the gravitino dark 
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matter with mass of order 0.1-1 TeV makes it possible to overcome such cosmological prob¬ 
lem [3]. For gravitino dark matter models, it was pointed out that the decay channels of 
gauge/Higgs bosons plus lepton dominate the total decay width under bi-linear RPV, which 
can produce observable high energy hadrons and leptons [4] . Another possibility is axino — a 
superpartner of axion, which in some supersymmetric scenarios becomes the LSP [5-10] and 
plays the role of dark matter [11-13] (see also Refs. [14, 15] for cosmology of supersymmetric 
axion model). Recently Ref. [16, 17] studied the axino dark matter with RPV in order to 
explain the observed baryon asymmetry. In this case RPV is crucial for baryogenesis, which 
is directly connected to hadronic products in axino decay. 1 Finally, anomalous cosmic ray 
positron fraction observed with PA ME T, A [19, 20], Fermi [21] and AMS-02 [22, 23] motivated 
intensive (in many cases phenomenological) discussions on leptophilic models, where dark 
matter particles decay (or annihilate) mainly into multiple leptons. 

In this paper, we study the impact of dark matter decays on the latest gamma-ray 
observations with the Large Area Telescope (LAT) onboard the Fermi satellite. In order 
to obtain stringent, yet robust, constraints, we adopt the latest data for the extragalactic 
gamma-ray background (EGRB) [24], which is isotropic gamma-ray radiation from every 
direction of the sky. Strictly speaking, EGRB is approximately isotropic since it partly 
consists of gamma rays from point sources (we will see in detail later). However, it is a good 
approximation to say that EGRB is isotropic statistically. Since we compare the EGRB data 
with the contribution from dark matter decays that is the statistical average of gamma-ray 
radiation from each of dark matter structures, this yields results that are largely independent 
of any uncertain inputs such as density profiles of dark matter halos. In addition, we take into 
account secondary gamma rays due to the inverse Compton (IC) scattering by the charged 
particles from decays off the CMB photons as discussed in Ref. [25]. 2 This again yields no 
theoretical uncertainties if a decaying dark matter model is specified; we do not include any 
other seed photons such as regular optical emission from galaxies. It is also known that 
hadronic decay of dark matter is strictly constrained by measurements of cosmic ray anti¬ 
protons [27-31]. Besides the fact that these anti-proton constraints are subject to much 
larger uncertainties due to Galactic cosmic ray propagation, we will show that the latest 
EGRB measurement gives more stringent constraints at this point. This should of course be 
revisited after the release of AMS-02 data for anti-protons and other light elements. 

Compared with the earlier work (see, e.g., Refs. [32, 33] and references therein) based 
on the Fermi-LAT EGRB data collected for 10 months [34], we improve upon the fol¬ 
lowing points. Firstly, we adopt the latest EGRB data measured with Fermi-LAT for 50 
months [24]. 3 Secondly, we address the systematic uncertainties that come from the mod¬ 
eling and subtraction of the Galactic foreground emission, and find that they reflect on 
uncertainties for lifetime constraints by a factor of a few. Lastly, in addition to purely phe¬ 
nomenological power-law modeling of other components, we use more realistic astrophysical 
components modeled individually with latest gamma-ray luminosity functions calibrated with 
other wavebands. We also adopt the uncertainties associated with these modeling as prior 
information in Bayesian statistical analyses to obtain lifetime constraints of dark matter. It 
is found that the obtained lower limit on the decay lifetime is improved by about an order 
of magnitude compared with Ref. [32], for all the hadronic and leptonic decay channels. For 
example, for a conventional model of decays into bb, we find that the lifetime shorter than 

1 This is an explicit realization of the baryogenesis, which is proposed in Ref. [18]. 

“IC photon in annihilating dark matter is studied in Ref. [26]. 

3 Constraint on annihilating dark matter in light of the 50-month data is studies in Refs. [35-37]. 
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10 28 s is excluded for dark matter masses between a few 100 GeV to 1 TeV. In addition, it is 
found that the parameter space which is motivated to explain the positron excess is almost 
excluded. 

This paper is organized as follows. In Sec. 2, we give formulation of the EGRB compu¬ 
tation from dark matter decays, and several models that give finite lifetime to dark matter. 
Section 3 then discusses three astrophysical components that are considered for the analysis. 
Results are presented in Sec. 4 along with details of the analysis, and we discuss several im¬ 
portant implications in Sec. 5. Throughout this paper, we use cosmological parameters given 
by Planck collaboration combined with WMAP [2]: present Hubble rate, and density divided 
by critical density of matter, dark matter, dark energy are Hq = 67.04 km s ” 1 Mpc ” 1 (i.e., 
h = 0.6704), = 0.3183, Odm = 0.2678, = 0.6817, respectively. The critical density is 

obtained by p c = 1.054 x 10 ” 5 h 2 GeV cm” 3 . 


2 Decaying dark matter 


Gamma rays from decaying dark matter is determined with three ingredients: dark matter 
mass, lifetime, and energy distribution of particles per decay, i.e., 

m dm , r dm , ^ (7 = 7 , e ± ,---) , ( 2 . 1 ) 

where the subscript ‘dm’ stands for dark matter. First two are free parameters, while the 
dNj/dE is computed for each dark matter model. We formulate the EGRB from decaying 
dark matter based on Refs. [4, 25], then discuss dark matter models on which we focus. 


2.1 Formulation 

Dark matter contribution to the EGRB consists of two components: primary and secondary 
gamma rays. Secondary gamma rays are produced by IC scattering off the CMB photons due 
to electrons and positrons from dark matter decay. The total gamma-ray intensity < h^ m (£' 7 ) 
is expressed as 

< 3 ^ m (£ 7 ) = [ dtQ^E^E'J , ( 2 . 2 ) 

47T m.d m 7"d m J 


where c is speed of light, t is cosmic time, E' = (1 + z)E 7 is the energy of 7 -rays when they 
are produced at redshift z, and 


Q^ m (£ 7 , E'J = e~ T ^\l + z) [Rprim(£ 7 ) + 7>ic (E^)] 


(2.3) 


where t(z,E 7 ) is optical depth, which we adopt the result given in Ref. [38]. V pr i m (Ely) and 
Vi c {Ej) are energy distributions of primary and IC 7 -rays, which are given by 


V 


prim 


Pic(7) = 


( 7 ) = ^(K 


dE 
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dE P dE. 


da\ C 


7CMB dE\ 


(K 1 E e , Eryg^ )/cmb(T 1 7cmb ) 


Y e (E e ) 
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(2.4) 

(2.5) 


respectively, where dcric/dE' is the differential cross section of the IC process, /cmb is the 
spectrum of the CMB (per unit volume per unit energy), and b lc = {A/‘i)a^{E e /m e ) 2 
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is the energy loss rate of the eb due to the IC process. Here gt is Thomson scattering cross 
section and Pqmb — 0.260 eVcm~ 3 , and Y e (E e ) is defined by using dN e ±/dE as 


fC 

Y e (E e )= j 
Je, 


dE 


dN e + , _ dN e - , _. 

(E) + —^(E) 


dE 


dE 


( 2 . 6 ) 


Then after changing variable from t to E' it is straightforward to obtain 


p2^dm/p \ c ^dm Pc Ey f°° , E y Q 1 (E^,E^) 

E ^ (E>) - L dE -K 
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(2.7) 


2.2 Dark matter models 


Decaying dark matter models can be classified into three categories by the products after 
cascade decay: (i) lepton, (ii) hadron + lepton and (iii) hadron. Decaying dark matter models 
of class (i) and (ii) are particularly popular since energetic positron due to the decay can 
account for the cosmic-ray positron excess which was reported by PAMELA [19, 20], Fermi- 
LAT [21] and AMS-02 [22, 23], with a proper choice of Td m and vn ^ To do comprehensive 
analysis, we consider several dark matter models in supersymmetric theory and compute 
gamma-ray distribution for each case. 

In supersymmetric models, the LSP becomes unstable once RPV operators are intro¬ 
duced in the superpotential, 



IML;H u + \ ijk QiLjb c k + b 3k L t L 3 El + A '^Dpi , 


( 2 . 8 ) 


where i,j,k are generation indices, Qi(Li) is left-handed quark (lepton), U E , D E (Ef) are 
right-handed up- and down-type quarks (lepton), and H u is up-type Higgs. We use ‘hat’ to 
represent chiral superfields (superpartners are expressed by using ‘tilde’ in the later discus¬ 
sion). The final state of dark matter decay depends on which of the above RPV interaction 
terms are operative. For the purpose mentioned above, we consider wino, sneutrino, gravitino 
and axino as the LSP. Simulated gamma-ray fluxes are summarized in Fig. 1, for which we 
will give the details below. 


2.2.1 Wino 

Winos are superpartner of W bosons. In a wide range of anomaly-mediated supersymmetry 
breaking scenario [39] , neutral wino W 0 is the LSP and it is a viable candidate for dark matter 
with a mass of a few hundred GeV [40] or around 3 TeV [41]. Besides, scalar partners of 
quarks/leptons are much heavier than TeV scale in anomaly mediation, which is consistent 
with 126 GeV Higgs boson discovered at the LHC [42, 43]. Testing wino dark matter is 
pointed out to be promising in direct [44] or indirect [45] detection experiments. More 
recently cosmic rays from decaying wino is studied in LLE C RPV [46] to explain the latest 
cosmic-ray positron excess reported in AMS-02 [22, 23]. We study the same type of RPV 
considered in Ref. [46] to simulate gamma ray. 
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Since squarks are much heavier than winos, it is instructive to write down interactions 
of W° which are relevant for the decay in higher dimension operators. In large slepton limit, 
they are given by 


E'mt — 


A' 


ijk 


A 2 


{W°P L ej)(e h P LVj ) + (W°P L ^)(e fe P L e,) 


+ h.c. 


(2.9) 


where P L = (1 - 7 5 )/2, 1/A 2 = y/2g 2 / , and summation of flavor index should be taken 

"7 

for i < j. g 2 is SU(2)l gauge coupling constant and is soft supersymmetry breaking 
mass for left-handed slepton. We have assumed rri^ < rn^ for i < j. Wino decay as 

W° —> ^ejef {piefe^) and Vje^ef [pjeje^). If lepton masses are ignored, then the energy 
distribution of final state leptons is simply given by 


d2 r W°-^e7e+ 

dzjdzk 

d 2 r W p^ Uje - e + 

dzidzk 


d 2 T l y 0 ^- i e+ e - 

dzjdzk 

d 2 T w 0_> gj . e + e - 
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5127T 3 A 4 


I A' 


ijk\ m w0 

5127T 3 A 4 


Zj(l - Zj), 

(2.10) 

Zj { 1 - Zj), 

(2.11) 


where m^o is wino mass and z a = 2 E a /m^ 0 (a = i. j, k) are defined by the energy E a of 
Vi, ej,ef an d so on - They satisfy Zi + Zj + Zk = 2 and 0 < Zij t k < 1 (see also Ref. [47]). Then 
the energy distribution for charged leptons in the final state is given by 


dN e . 

dz e . 


9 / \ 

12z 2 (l-z e .), 


dz. 


— 2 z e k (3 2z efc ) , 




( 2 . 12 ) 


in a single process W° -» vie - ef, and 


dN e(j 

dz ep 


2^e^(3 ^ z ep) i 


(2.13) 


where /3 = i,k in IT 0 —> Vje^ef. Energy distribution for each particle in CP conjugated final 
state is the same as in original state. For example, if A122 is only relevant, then the final state 
in wino decay is v e g~fi + (v e g + g~) and v^e~g + (v ll e + g~). Thus no primary gamma rays or 
(anti-)protons are generated. In such case, IC photons from eP (which includes electron and 
positron from /r^) are the only observable gamma rays. If is produced, then its cascade 
decay produce energetic gamma rays. We use PYTHIA [48] to compute dNj/dE (/ = 7, e ± ) 
in the primary decay of dark matter. 


2.2.2 Sneutrino 

As a reference, we also consider right-handed sneutrino dark matter in LLE C RPV. Supposing 
that neutrinos are purely-Dirac fermions, neutrino Yukawa couplings are very small, i.e. 
C>(10 -13 ). Since right-handed sneutrinos, superpartners of right-handed neutrinos, interact 
with the other particles via the Yukawa couplings, the lightest one can be dark matter if it 
is the LSP [49]. In LLE C type RPV, right-handed sneutrino vm decays to charged leptons 
if Ik - We will compute gamma-ray flux for the final state Z^Z/"- 4 

4 Electroweak correction to the final state becomes important when dark matter mass is extremely large. 
However, we have checked numerically that it does not affect the result when the mass is less than lOTeV. 
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2.2.3 Gravitino 


In imparity conserved case, gravitino LSP is typically disfavored in cosmology when it has 
electroweak-scale mass. This is because the next-LSP (NLSP), which is usually the LSP 
in the minimal supersymmetric standard model and decays to gravitino and standard-model 
particles via Planck-suppressed interaction, becomes long-lived and may decay after big-bang 
nucleosynthesis (BBN) started. In the RPV case, however, NLSP can decay to standard- 
model particles via RPV interaction before BBN starts. Thus NLSP decay does not affect 
BBN. Furthermore, gravitino is cosmologically long-lived since the decay rate is suppressed 
by the Planck mass and a small violation of R-parity, thus it can play the role of dark 
matter [3]. In bi-linear RPV, for example, gravitino t/ 3/2 decays to 7 z/j, W±lf, Zvi and 
hut if kinematically allowed. It is shown that the branching fraction of ^ 3/2 Iis the 
largest [4], which means that the decay products after cascading decay of W±lf are mixture of 
high-energy electrons/positrons, gamma rays and protons/anti-protons. 5 Recently decaying 
gravitino dark matter is revised under bi-linear RPV [50, 51] to account for cosmic-ray 
positron excess observed in AMS-02 experiments. In our analysis we simply consider ip ^/2 —> 
W±lf to avoid introducing parameters involved in determination of each branching fraction, 
and simulate primary and IC gamma rays. Inclusion of the other channels is straightforward 
(see Ref. [4]), and that will not change our result significantly. 


2.2.4 Axino 

Axino, a superpartner of axion originally introduced to solve the strong CP problem [52], 
can be dark matter in some supersymmetry models. In addition to that, if baryon number 
is broken by U C D C D C term with \"- k < l, 6 baryon asymmetry is generated by moduli decay 
before BBN to give the observed baryon number [17]. Even in such a large RPV, it is shown 
that axino with a mass of 0(10 GeV) can be long-lived to be dark matter [16]. 

In this model, scalar superpartners are much heavier than the electroweak scale. The 
interactions of axino (denoted as a) with squarks are then written in higher dimension oper¬ 
ators. Ignoring down-type squarks and left-right mixing in squark sector for simplicity, it is 
given by 

X" 

Ant = -^(uiP L a)(d k P L d%) + h.c. , (2.14) 


where 1/A 2 = 2 g^s/rri-. . Here m,u R . is supersymmetry breaking mass of ur ?: , and is 
one of the dimensionless couplings of axino-quark-squark interaction, which is suppressed 
by axino decay constant (see Refs. [12, 55]). Color indices are implicit. Then axino decays 
in three-body process: a —> Uidjdk, Uidjdk■ If quark masses are ignored, then the energy 
distribution of final state quarks is simply given by 

d 2 Ta^ Ui djd k d ^a—>Uidjd k 3|A i fc | m| , , 

dzdZi = ifeitfej = l 28 ^A^ i(1 - a) - (2 ' 15) 

where is axino mass, Zi = 2Ei/rria and Zj = 2Ej/m.a are defined by the energy Ej and Ej 
of m and dj, respectively. Then it is straightforward to give the energy distribution of each 

5 Neutrinos are also produced, but they are irrelevant in our discussion. 

6 See Refs. [53, 54] for realization of such case. 
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E 7 [GeV] 


Figure 1. Gamma-ray fluxes from various decaying dark matter (rridm = 1 TeV, 7d m = 3 x 10 27 s). 
Plots give fluxes from decay channels: (a) (i/ ey u + /r _ ) and fi + (b) /i + fj ,~, (c) 

t + t~, (d) W-/E, (e) uds ( uds), and (f) bb. Data points with error bar and a band of the EGRB 
observed by Fermi-LAT is also shown [24] (see Sec. 3). 


final sate quark: 

9 . . d.N q 9 . . . 

^ = 12 ^( 1 -*), -^ = 2 %$- 2 Zj ) , ( 2 . 16 ) 

in a single process a —>• Uidjdk■ The energy distribution for dk is the same as dj. These quarks 
are hadronized to produce mesons, which decay to gamma rays and electrons/positrons, 
and electrons/positrons become source of IC photons. In later numerical analysis, we also 
compute a case of final state bb for comparison, which would be useful for those who are 
interested in. 


2.3 Gamma-ray fluxes (examples) 

In Fig. 1 gamma-ray fluxes in various decaying dark matter models are plotted. For lep- 
tophilic case, result is shown for a case where only A 7 122 is relevant (dubbed as “zxZ + Z _ ”) in 
W° dark matter, while decay channels and t + t~ are considered in vr decay. It is seen 

that the gamma-ray spectra from LLE C and /i + ^ are quite similar. On the other hand, in 
t + t~ , the spectrum has double peaks. This is due to primary gamma rays produced from 
cascade decay of tau, which gives another gamma-ray flux in high energy region. For hadron- 
ically decaying dark matter, the axino decay via A'/ 22 is considered (denoted as “uds”). The 
spectrum shows similar behavior to t + t~ case and bb channel as well. Finally, the flux from 
decaying gravitino to W ± /x =F is expected to have a property in the middle of leptophilic and 
hadrophilic cases, which is in fact seen in the figure. 
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Figure 2. EGRB spectrum from blazars (solid), star-forming galaxies (dashed), misaligned active 
galaxies (dot-dashed), and dark matter decay (dotted; ma m = 1 TeV, Td m = 3 x 10 27 s, bb decay 
channel). Data with error bars show the latest Fermi-LAT measurement [24], and a band attached 
to them show systematic uncertainty due to foreground subtraction (see Sec. 4.2). Note these data 
are total EGRB, which includes already resolved sources. Unresolved EGRB (or also referred to as 
isotropic gamma-ray background), on the other hand, is about a factor of two smaller [24]. 


3 Astrophysical source 

Conventional astrophysical sources emit gamma rays through particle acceleration followed 
by interactions between accelerated charged particles and surrounding media or photon fields. 
The total EGRB should be made of both dark matter and astrophysical components: 

4> 7 (F 7 ) = <^ stro (F 7 ) + ^ m (F 7 ) . (3.1) 

In this section we discuss three possible astrophysical sources in extragalactic region. They 
are blazars, star-forming galaxies (SFGs), and misaligned active galactic nuclei (rnAGNs). 

3.1 Blazars 

In the GeV gamma-ray sky, the most dominant source is blazars, which are one class of 
active galaxies whose jets are directed towards us. More than thousand blazars have been 
identified with Fermi-LAT [56] and they indeed make up a significant fraction (~70%) of the 
total measured EGRB intensity [24, 35]. Therefore, it is essentially important to include this 
source class with a particular care. 

Since many blazars have been identified with known redshifts, it is possible to recon¬ 
struct its gamma-ray luminosity function (i.e., number density of blazars per unit luminosity 
range). This has been done extensively, and thus far, the luminosity-dependent density 
evolution model that is motivated by possible underlying correlations between X-ray and 
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gamma-ray luminosities gives the most reasonable fit to both the luminosity and redshift 
distributions. The luminosity functions constructed this way can be found in Refs. [57, 58] 
for flat-spectrum radio quasars and BL Lactae objects, respectively (each is a subclass of the 
blazar population). In this study, we adopt the EGRB spectrum due to blazar population 
from Ref. [35] based on the latest reconstruction of the blazar luminosity function, which 
is shown in Fig. 2. The current uncertainty range of the spectrum is about ~30%, almost 
independent of energy. Fig. 2 shows that the blazars makes up of significant fraction of the 
EGRB. This is not very surprising because the EGRB data shown in the figure also includes 
all the resolved extragalactic sources, most of which are blazars. However, it is quite clear 
that they cannot be the only source that can explain the EGRB in the entire energy region 
between 100 MeV and 820 GeV. 

3.2 Star-forming galaxies 

SFGs are galaxies that are undergoing active star formation just like our own Milky-Way 
Galaxy. Even though they are in general much dimmer than the blazars, already several of 
them are detected with Fermi-LAT [59]. The same study also found a tight correlation be¬ 
tween the gamma-ray luminosity and infrared luminosity for these identified galaxies. Assum¬ 
ing that this relation universally holds and by using the precise infrared luminosity function 
of SFGs, one is able to construct the gamma-ray luminosity function of SFGs. Reference [60] 
computed the EGRB intensity due to SFGs, starburst galaxies, and SFGs containing AGNs 
separately, by adopting the latest infrared luminosity function constructed with the Herschel 
telescope [61]. The same paper found that the uncertainty of the EGRB spectrum came 
mainly from those due to gamma-ray-infrared luminosity correlation, and it yielded about 
~60% uncertainties. It also yielded larger flux than the previous model [59]. Fig. 2 shows 
the SFG contribution (that also includes contributions from starbursts and SFGs containing 
AGNs). The peak around ~0.3-0.5 GeV comes from kinematical argument of the pion de¬ 
cay that is produced in the interactions between cosmic rays and interstellar medium. The 
spectrum becomes suppressed above ~20 GeV because of the effect of gamma-ray absorption 
due to the extragalactic background light. The SFG spectrum nicely complements that of 
blazars especially at low energies. 

3.3 Misaligned active galactic nuclei 

Yet another population, although much less certain, is the mAGNs, whose jets are not di¬ 
rected towards the Earth. (They are also known as radio galaxies.) Even though they are 
less bright than the blazars in gamma rays, a few of them are identified with Fermi-LAT. The 
gamma-ray luminosity function for this population is constructed with the radio luminosity 
function and measured correlation between radio and gamma-ray luminosities. Although 
with much larger uncertainties, this has been done in the literature [62, 63], and for this 
work, we adopt the model of Ref. [62] that is also shown in Fig. 2. The uncertainty is still 
quite large due to paucity of the detected sources, which is typically a factor of three. 

4 Results 

We give the lower bound for the lifetime of dark matter. Based on the discussion in Sec. 2, we 
consider following final state in dark matter decay; u e ^x~fi + (v e n + n~) & n^e~fi + (u fl e + n~), 
t + t~, W ± n T , uds ( uds ), bb. 
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Table 1 . Priors for dark matter and astrophysical parameters. 


Prior 

7dm 

/blazar 

/sFG 

/mAGN 

Normal 

20 < log(r dm /s) < 30 

'Tog / = 0-1 

‘Tog/ = 0'“ 

‘Tog/ = 0-5 

Flat 

20 < log(r dm /s) < 30 

-5 < log/<2 

-5 < log / < 2 

-5 < log / < 2 


4.1 Analysis 

With the EGRB data measured with Fermi-LAT [24] and both the dark matter and as¬ 
trophysical models discussed in the previous sections, we perform the Bayesian statistical 
analysis using the Markov-Chain Monte Carlo simulation to obtain 95% credible lower limit 
on dark matter decay lifetime for given dark matter model characterized by mass and decay 
channel. 

The posterior probability distribution of parameters d given data d , P(d\d), is related 
to the prior P(d) and the likelihood function C(d\d) through 

P(d\d) oc P(d)£(d\d) . (4.1) 

We assume that the likelihood function £(d\d) is approximated to be normal distribution 

✓V^-TT 1 { [^-^7 m (^kdm)-$f r °(^|/blazar,/sFG,/ m AGN)] 2 

11 ’ B V2itr t< eXP ( 24, 

(4.2) 

where Ej and 4>j represent the energy and the EGRB intensity data in i-th bin, respectively. 
The lcr errors for 4>j, assumed to be uncorrelated between different energy bins, are shown 
as (j$ ]; (and also as the error bars in Fig. 2). For theoretical parameters, we adopt the dark 
matter lifetime, Td m , and energy-independent normalization for each of the astrophysical 
components, /blazar j /sfg> and / m AGN> for which / = 1 recovers our reference spectra shown 
in Fig. 2. 

For the prior P( th m ). we adopt a flat distribution in logTd m in the range of 20 < 
log(rd m /s) < 30, which covers quite a wide range without any bias. Regarding the prior 
for the astrophysical parameters, we adopt normal distributions centered on log / = 0, and 
with standard deviation of <Rog/ blazar = 0.1, o'log/sFG = 0.2, and 01 O g/ mAGN = 0.5, based on 
the errors 25%, 60%, and 200%, for the blazars, SFGs, and mAGNs, respectively. Although 
the errors of blazar depend on the energies of gamma ray, they are ~25% in most energy 
bins, while they are larger (~30%) at a few high-energy bins. These ‘Normal’ priors, also 
summarized Table 1, accommodate the best knowledge thus far obtained for the astrophysical 
sources as discussed in Sec. 3. We note that since the information of the resolved sources are 
already used in order to obtain these Normal priors, here we adopt the errors on unresolved 
EGRB intensity [24] in denominator of Eq. (4.2). For the gamma-ray intensity, on the other 
hand, we use the data of the total EGRB. 

In order to give more conservative limits, we also discuss the case without any astro- 
physical prior information, for which we simply adopt the log-flat priors for all / parameters, 
between —5 < log / < 2. These ‘Flat’ priors are summarized in Table 1. Since / parameters 
follow the log-flat distribution down to a very small value (/ = 1CF 5 ), the parameter space 
includes a case where the astrophysical sources give virtually zero contribution. Thus, in 
order to analyze consistently, we adopt the data and errors a$ i for the unresolved EGRB 
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Table 2. Priors for phenomenological parameters for ‘astrophysical’ component. 


Prior 

/pl 

7 

Ecut 

PL 

-2 < log /<2 

-3.5 < 7 < -1 

100 < £ cut /GeV < 600 


(see caption of Fig. 2). We note that this approach is somewhat pessimistic, because many 
astrophysical sources are known to give contributions to the unresolved part of the EGRB, 
and hence taking extremely small values for / will result in significant underestimate of these 
astrophysical components. 

For purely phenomenological purpose and also in comparison to the results with 10- 
month Fermi data in the literature [32], we also adopt a single cutoff power-law component 
for non dark matter EGRB. Although from Fig. 2 such a single cutoff power-law model 
appears to give reasonable fit to the Fermi data (indeed it does [24]), it is no longer strongly 
motivated since the resolved components (mainly blazars) are already shown to feature quite 
different spectral shape than the EGRB data. For this phenomenological power-law model 
(labeled as ‘PL’), we adopt the normalization factor /pl, where /pl = 1 corresponds to the 
best fit parameter to the data on its own, the spectral index 7 , and the exponential cutoff 
energy E cut , following Ref. [24], For /pl, we use the log-flat prior, and for the other two flat 
priors in linear scale. The ranges of these priors are summarized in Table 2. 

4.2 Constraints on dark matter lifetime 

Figure 3 shows 95% credible lower limits on the decay lifetime as function of its mass for 
various decay channels. Here we adopted the astrophysical background models with Normal 
prior (Table 1), which should be regarded as our canonical results based on our most up- 
to-date knowledge on astrophysical sources. Thick solid curves in Fig. 3 shows the lifetime 
constraints obtained with the data shown in Fig. 2. These data are obtained after the 
Galactic foreground emission being subtracted, after it was modeled with a parameter set of 
the Galactic cosmic ray propagation (model A in Ref. [24].) We find that the new lifetime 
constraints are much more stringent than those obtained with the 10 -month data by up to 
more than one order of magnitude. 

In leptonically decaying dark matter models, such as and v e [i ~& 

/i + constraint on the lifetime becomes stringent almost monotonically as 

mdm and rdm < 10 27 -10 28 s around ru-dm ~ 10 TeV is excluded. This can be understood 
since gamma rays in IC process is the only contribution from dark matter and the peak of 
the intensity (i.e. E 2 ^ 111 ) shifts to the region where Fermi-LAT data have better accuracy 
as ?7idm increases. (See Fig. 1.) In the other cases, where decay products contain lots of 
hadrons, the lifetime shorter than 10 28 s is almost excluded for dark matter masses between 
a few 100 GeV and ~1 TeV. As opposed to the leptophilic case, the constraints get weaker for 
?7idm > 1 TeV. This is because the gamma-ray spectrum have two peaks which comes from 
the cascade decay and one from IC process. When 777d m > TeV, a “valley” in the spectrum 
enters in the range 1 GeV < E^ < 100 GeV, where Fermi-LAT data is given with small 
errors. Thus the constraint gets relaxed.' 

Quantitative arguments depend on the foreground model adopted and subtracted from 
the total gamma-ray emission. To this end, we repeated the same computation by using two 

7 Those interpretation is qualitative. Quantitative behavior is determined by non-trivial interplay between 
modeling of the foreground subtraction and choices of the prior. See discussion below. 
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Dark matter mass m dm [GeV] 


Figure 3. 95% credible lower limits on dark matter lifetime Td m as function of dark matter mass m d m , 
for decay channels: (a) v e n~ /i + (P e /r + jU _ ) and n + (b) (c) t + t~ , (d) W ± ^ T , (e) 

uds (uds), (f) 66. Astrophysical background models with Normal priors are adopted (Table 1). Thick 
solid, dashed, and dotted curves correspond to the EGRB data with different foreground modeling 
discussed in Ref. [24] (their models A, B, and C, respectively). Thin solid curve shows the lower limits 
obtained with the 10-month Fermi-LAT data [34] and the phenomenological power-law background 
modeling. 


different foreground models, B and C adopted also in Ref. [24]. Models A—C nicely covers 
regions shown as uncertainty band in Fig. 2. The dashed and dotted curves are the results 
corresponding to models B and C, respectively. This shows that the foreground modelings 
give uncertainty on lifetime constraints by about a factor of a few. 

The results of more conservative approach with Flat priors in Table 1 are shown in 
Fig. 4. As expected, in most cases, they are weaker than the ones with Normal priors (as 
shown in Fig. 3) by about a factor of a few. Exceptions are at high dark matter masses 
for (c)-(f), where they give stronger constraints; this is likely caused by interplay between 
different choices of priors and the data (the total EGRB data for the Normal priors, while 
the unresolved EGRB data for the Flat priors). 

In order to compare our results with the previous ones in the literature (e.g., Ref. [32]), 
we also computed the lifetime constraints by using the 10-month Fermi-LAT data [34], Here 
we modeled the other background component as a single power law (Table 2), and the re¬ 
sults are shown as a thin curve in each panel of Figs. 3,4 and 5 for reference. Although the 
statistics adopted here is different than that in Ref. [32] (Beyesian versus frequentist), our 
results are in good agreement with theirs, proving the consistency of both the approaches. 8 

s The result for t + t~ in high mass region is different from Ref. [32]. This is because they used both 
published and preliminary data for F 7 > 100 GeV (at that time) while we use the published 10-month data 
only. In t + t~ case, gamma-ray spectrum from cascade decay is hard and the peak of the intensity is out of 
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Figure 4. The same as Fig. 3, but for the astrophysical background models with Flat priors (Table 1). 
Note that these priors are very conservative. 
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Figure 5. The same as Fig. 3, but for phenomenological background models with PL prior (Table 2). 


data region when mj m > TeV. Then the constraint becomes weaker. 
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Figure 6. The same as Fig. 3, but taking into account only dark matter contribution. 


In Fig. 5, we show results by using the phenomenological cutoff power-law model as the astro- 
physical component. Again, this is just for reference purpose, since such a single-component 
astrophysical modeling is no longer valid. 

Finally in Fig. 6 we repeated our analysis by switching off all astrophysical components, 
i.e. taking into account only dark matter component, for reference. This result corresponds 
to the most conservative limit on the lifetime, which may be helpful for some readers. 

4.3 Implications for particle physics models 

As discussed, dark matter models in which large amount of high-energy electrons and positrons 
are produced are phenomenologically motivated in order to explain the observed cosmic-ray 
positron excess. Typically a lifetime of 10 26 s with a mass of ~1 TeV is required to account for 
the excess. Our result indicates that such dark matter models are almost excluded already. 

The parameter region of decaying gravitino via LH U RPV, 7713/2 = 1—2 TeV and r 3 / 2 — 
10 26 s [51] to explain the latest AMS-02 positron data, is excluded in both Normal and Flat 
priors (Figs. 3 and 4). In Ref. [51] constraint from the EGRB is also discussed. Taking 
into account the astrophysical sources and 50-month Fermi-LAT data, they concluded that 
the decaying gravitino with the positron excess-motivated region is not excluded, which is 
inconsistent with our result. 9 The decay channel fj, + fi~ to explain the positron excess is 
excluded as well. The same conclusion is already stated by Ref. [32]. We have confirmed 
their conclusion in more sophisticated treatment for astrophysical sources. Lastly decaying 

9 This statement was based on the first version of Ref. [51]. In the second version they fixed a bug in the 
computation of gamma-ray flux and gave a consistent result with ours. We thank German Gomez-Vargas for 
pointing this out. 
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wino dark matter via LLE C interaction to give sufficient positron flux to explain AMS-02 
data is also excluded for thermal wino with a mass of 3 TeV. 

The lifetime of dark matter which decays mainly through U C D C D C RPV is severely 
constrained too. In such case, though decaying dark matter no longer explains the positron 
excess, decaying axino, for instance, with a mass of 0(10 GeV) is motivated in the context 
of baryogenesis, which has been mentioned in the previous sections. The results obtained in 
this study shows that = 10 26 -10 27 s for 0(10 GeV) axino mass are partly excluded. This 
means that the parameter space of axino dark matter in moduli-induced baryogenesis may 
be probed in the future observation. 

Here we have discussed the implication of our numerical results for the dark matter 
models described in Sec. 2.2. However, it is obvious that the results given in Sec. 4.2 does 
not depend on what dark matter is. Only the decay products of dark matter, i.e. the final 
state, change the constraints on the lifetime of dark matter (as function of its mass). Thus, 
if gravitino decays via LLE C RPV, for example, then the result of (a) in Figs. 3 and 4 can 
be roughly applied. 10 


5 Discussion 


5.1 Comparison with constraints from galaxy clusters and dwarf galaxies 

There are other regions of the sky that can be used to constrain dark matter decay through 
gamma rays. The largest virialized dark matter structures, known as clusters of galaxies, 
were studied in the literature for the purpose [64, 65]. Dwarf spheroidal galaxies are nearby 
dark matter substructure in the Galaxy that are highly dark matter dominated, and hence 
they are considered ideal for dark matter searches. The dwarf galaxies were also studied for 
constraining dark matter decay [64] . Here, we compare them with the EGRB that was used 
in our study. 

Gamma-ray intensity <h 7 is proportional to line-of-sight integral of dark matter density 
as was discussed in Sec. 2.1. In case of the EGRB where the integration should be taken up 
to the Hubble horizon, this quantity is 


•Irgrb 


cDdm Pc f ^ 

Ho J 4" f^m(l + z)3 


(5.1) 


and is 3 x 10 22 GeV cm -2 if we integrate over 0 < z < 5. Note that we here ignored gamma- 
ray absorption as well as the gamma-ray spectrum per decay dN^/dE found in Eq. (2.7), on 
order to have order-of-magnitude insight of the relevant term. This is now compared with the 
same quantity for the galaxy clusters and dwarf galaxies, after smeared over the Fermi-LAT 
angular resolution. Reference [64] computed it for both the sources, and found that they 
were typically 2 x 10 22 GeV cm -2 for the most promising cluster (Fornax) and more than one 
order of magnitude smaller for all the known dwarf galaxies, after averaging over the region 
with 1° radius. They are not only smaller than Jegrb but also subject to uncertainties on 
dark matter density profiles. In addition, the EGRB spectrum is inferred from the all-sky 
data, where there are many more photons available for the analysis than small regions around 
galaxy clusters or dwarf galaxies. Therefore, we conclude that the EGRB provides the best 
opportunity to constrain dark matter decay with gamma rays. 

10 Strictly speaking, the energy distribution of charged leptons from gravitino decaying in three body are 
different. However, the difference will not change the constraints on the lifetime significantly. 
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5.2 Comparison with constraints from Galactic anti-protons 

If dark matter decays hadronically, it produces lots of protons and anti-protons. Since high- 
energy anti-protons are rare in canonical Galactic activities, measurements of cosmic-ray 
anti-proton constrain hadronic decay of dark matter. However, it is known that there is huge 
theoretical uncertainty in the calculation of cosmic-ray anti-proton in the Galactic region. 
Cosmic-ray anti-proton flux can be calculated in a Galactic propagation model, which has 
several parameters. They are determined by using the measurements of other cosmic rays, 
such as proton flux and boron-to-carbon ratio, etc., which unfortunately cannot determine 
all of the parameters. Propagation models called ‘MAX’, ‘MED’ and ‘MIN’ are the typical 
examples [66]. Although they are consistent with the measurements of other light elements 
in cosmic rays, they give different predictions for anti-proton flux. In decaying dark matter 
model, the anti-proton flux changes up to an order of magnitude between MAX and MIN 
models (see Ref. [28], for example). 11 In addition, the constraint becomes much weaker in 
low dark matter mass region, i.e. md m 100 GeV [28]. On the contrary, the constraints 
from the EGRB we have studied in this paper have less uncertainty compared to those in 
cosmic ray p. It is noticed that there also exists an uncertainty which is involved in Galactic 
propagation model. The 50-month data from Fermi-LAT [24] is given by using GALPROP 
package [69] in order to simulate diffuse gamma rays in the Galactic region. They discuss 
the uncertainty by computing inner-Galactic gamma rays in different Galactic models, which 
corresponds to models A, B, and C described in the previous sections. Our results show 
that the constraints on dark matter model changes within an order of magnitude in the 
three models. Thus observation of the EGRB gives much more robust limits on decaying 
dark matter model with less uncertainty compared to the limits from anti-proton flux. In 
addition, it should be noted that the 50-month data gives much more stringent constraint on 
the lifetime than 10-month data in all types of dark matter models studied here. Thus more 
accumulated data with expanded analyses including anisotropies [70, 71] will give probably 
the best limits on decaying dark matter models in the future. 

6 Conclusions 

We have studied the extragalactic gamma-ray background in decaying dark matter models. 
Decaying dark matter produces gamma rays both primarily from the cascade decay and 
from electron/positron through inverse Compton scattering off the microwave-background 
photons. If the final state, along with its mass and lifetime of decaying dark matter, is 
specified, the gamma-ray spectrum can be computed with little theoretical uncertainty. As- 
trophysical sources of gamma rays in extragalactic region, on the other hand, have been 
resolved with gamma-ray telescopes such as Fermi. Main components, which are blazars, 
star-forming galaxies and misaligned active galactic nuclei, give rise to consistent flux with 
the measurement of the gamma-ray background according to Fermi. We have used the up¬ 
dated astrophysical models for those components and the latest Fermi data after 50-month 
sky survey to constrain decaying dark matter models. In order to make qualitative statement 
for various class of models, we have considered supersymmetric dark matter in R-parity vi¬ 
olation. Supersymmetric model accommodates variety of dark matter candidates, such as 
wino, gravitino, sneutrino and axino and so on. They decay to standard-model particles 
leptonically or hadronically, depending on the relevant R-parity violating interaction. We 

11 See also Refs. [67, 68] for recent development. 
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have calculated gamma-ray flux in the models where dark matter mainly decays to l~j~l k , 
Vieje k (vi e+e- k ), WHf, Uidjdk ( Uidjdk ), and q^. The numerical analysis shows that the 
constrains on the lifetime of dark matter becomes more stringent by a factor of a few to an 
order of magnitude, compared to the past work. As the result, we found that leptonically 
(and hadronically) decaying dark matter which explains the cosmic ray positron excess is 
excluded in most cases. Hadronically decaying dark matter is also severely constrained. For 
example, the current Fermi data are beginning to exclude part of parameter space of decaying 
axino dark matter in moduli-induced baryogenesis, which indicates that future observation 
of the extragalactic gamma ray may probe this model. 
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